Collective effects of multi-scatterer on coherent propagation of photon in a two 

dimensional network 
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We study the collective phenomenon in the scattering of a single-photon by one or two layers of 
two-level atoms. By modeling the photon dispersion with a two-dimensional (2D) coupled cavity 
array, we analytically derive the scattering probability of a single-photon. It is discovered that in the 
case with one layer of atoms, the atomic collective excitation leads to a shift for the single-photon 
scattering spectrum. Such a shift is related to the density of the atomic ensemble. For the case 
with two layers of atoms, an inter-layer effective coupling appears and induces an electromagnetic- 
induced-transparency-like phenomenon for the single-photon scattering. Our result provides a new 
scheme of analyzing photon coherent transport in 2D and may help to understand the current 
experiments about the high energy photon scattering by the layer nuclei material. 

PACS numbers: 32.70.Jz, 32.80.Qk, 42.50.Nn 



I. INTRODUCTION 

For the purpose of controlling the transport and the 
scattering of a single-photon in a full quantum fash- 
ion, progress has been made in confined one-dimension 
systems[l-3]. The transport of a single-photon can be 
coherently manipulated by the interaction between the 
photon and the doped natural or artificial atoms. The 
corresponding physical implementations could be real- 
ized in several different ways, such as: defected pho- 
tonic crystals i4ft8i|, superconducting transmission line 
resonators (9l-[l3| and coupled cavity arrays (CCA) [Tsl - 
\it\ . By modeling the photon hopping between the cav- 
ities as a tight-binding Hamiltonian, the photon in the 
cavity array displays a nonlinear dispersion relation. This 
fact makes the cavity array a well-controllable platform 
for studying the photon transport with different quantum 
natures. For example, with a tunable two-level atom in- 
side one of the cavities, the transmission and reflection of 
a single-photon can be well controlled in one-dimensional 
CCA 

In order to fabricate integrated all-optical on-chip 
devices, it is necessary to move forward to the two- 
dimensional (2D) structures. Based on CCA systems, 
a lot of theoretical investigations have been proposed to 
use this system as an alternative for quantum simulation 
of many body physics (T8l - [2^ . However, the studies of 
coherent scattering of a single-photon in this 2D cavity 
array is still lack to our best knowledge. 

In this paper, we generally study the scattering prob- 
lem of a single-photon in the 2D cavity array with two- 
level atoms. In our model, the atoms are periodically 
located in one or two rows of the cavities. We analyti- 
cally solve the scattering problem, and obtain the exact 



expression for the single-photon scattering probability. 
These results explicitly display the collective coherent ef- 
fects in the scattering process. For the case with one 
layer of atoms, we find that the photon-atom interac- 
tion induces a collective energy shift for the atom, and 
the scattering probability takes the maximum value when 
the incident photon is resonant with the shifted atomic 
energy. In the presence of two layers of atoms, we find 
that the effective couplings between atoms in different 
layers are induced by the photon-atom interaction. As 
a result, the atoms are in two dressed states with dif- 
ferent effective energies, and the maximum of the scat- 
tering probability appears when the incident photon is 
resonant with one of these two states. Similar as the 
atomic susceptibility in the system with electromagnetic 
induced transparency (KIT), the scattering probability 
has a double-peak behavior as a function of photon-atom 
detuning. Therefore, the single-photon scattering in 2D- 
CCA implies more abundant phenomenon than in the ID 
case, and this simple model may offer us an intuitional 
microscopic understanding of the recent experiments of 
X-ray quantum optics [2^, 2M ■ 



'Electronic address: |pengzhang@ruc.edu.cn| 



The rest of the paper is organized as follows. In Sec. II 
we investigate the single-photon scattering in a 2D-CCA 
with one layer of atoms, and discuss the collective shift of 
the atomic energy. In Sec. Ill we consider the case with 
two layers of atoms, and illustrate the EIT-like behavior 
of the scattering probability. There are some conclusion 
and discussion in Sec. IV. Some detail of our calculation 
are presented in the appendix. 



2 



II. SINGLE-PHOTON SCATTERING WITH 
ONE LAYER OF ATOMS 

A. System and Hamiltonian 

We consider a 2D array of identical single-mode cav- 
ities as shown in Fig. [TJ We further assume that the 
photons can hop between neighbor cavities. Then the 
cavity array is described by a tight-binding model as 



+ 00 

X,y — — OG 



(1) 



Here ujc is the frequency of the photons in the cavities, 
^ is the hopping intensity or the inter-cavity coupling 
strength, a(a;,y) and a^^ are the annihilation and cre- 
ation operators of the photon in the cavity at position 
{x,y), respectively. Here and after we set h = 1. When 
there is a single photon propagating in the system, the 
eigenstate \k) of He takes the form of 2D plane wave with 
momentum k — [k^, ky): 
+00 



\k) 



1 

2^ 



E 

x,y — — oo 



vac) 



(2) 



with |vac) the vacuum state of all the cavities. The single- 
photon dispersive relation of the 2D cavity array 

H = '^(k:.,ky) = - 2$(cos k^ + cos ky) (3) 

is naturally obtained by the stationary Schrodinger equa- 
tion Hc\k) — e^\k). 

To explore the scattering character of the incident pho- 
ton on a collection of identical atoms with geometrical 
configuration, we embed a two-level atom in every d cav- 
ity along the y-axis, and the j-th atom is in the cavity at 
(0, dj). The Hamiltonian of these atoms is 



(4) 



with Ua the energy-level spacing of the two-level atom 
and \e)j {\g)j) the excited (ground) state of the j-th 
atom. The atom-photon coupling is described by the 
Jaynes-Cummings Hamiltonian 



V = ^(f7a(o,dj) |e)j- {g\ + h.c), 



(5) 



with ri the coupling strength. 

The setup studied in this paper can be physically im- 
plemented by placing the two-level atoms in the defected 
cavities of the 2D optical crystals [3. At the match 
condition (fca;(j/) '''/2)) the energy eg becomes a linear 
function of the momentum k as 6g ~ lOq -f- l^l^kx ky) . 
Therefore, in that region our model can also character- 
ize the scattering of photons on the array of atoms in 
the free space, e.g., the scattering process in the recent 
nuclear quantum optical experiments with X-ray. 




Figure 1: (Color online) The single-photon scattering by one 
layer of atoms in a 2D CCA. The two-level atoms are confined 
in one row of the array with periodic d (in the figure we show 
the case with d — 3). During the scattering process, the inci- 
dent photon (solid line) can be scattered in several different 
outgoing directions (waved lines). 



B. single-photon scattering state and T- matrix 

Now we calculate the scattering probability of a single- 
photon scattered by the two-level atoms in our system. 
To this end, we first derive the single-photon scattering 
state and the on-shell element of the T-matrix in this 
subsection. With the help of these results, we will ob- 
tain the single-photon scattering probability in the next 
subsection. 

The scattering state is given by the Lippman- 

Schwinger equation 



|vI/(+)) = |fc)|.g) 



{Ha + He) + iO+ 



V\¥+^) (6) 



with \g) = Ylj \9)j the collective ground state of all 
atoms. It is clear that in our system the total excitation 
E^r=-ooa(x,y)a(^,a) +Ei |e)i (e| is conserved. Thus, in 
the subspace with one excitation we can expand the sta- 
tionary eigenstate as 



= 1-^)1.9) +E/3.|vac)|g,). 



(7) 



Here \ip) is a single-photon state of the cavity modes, and 
e)^ is defined as \ej) = \e)j ® Hz^itj Iff); which represents 
the state with only the j-th atom excited. It follows from 
Eq.® and Eq.0 that |0) and f3j satisfy 



|vac) 



(vac I a 



(8) 
(9) 



Eqs. ([HI [H]) can be analytically solved by the follow- 
ing approach. First, we notice that our system is in- 
variant under the translation for d cavities along the 
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y-axis. Namely, the total Hamiltonian Ha + He + V 
of our system is commutative with the translation op- 
erator _D, which is defined as Da^^^y-jD'' = 0'(x,y+d)-, 
D |e,-) {e,\D^ = |e>i) (e,+i|, and D \~g) {~g\D^ = \~g} {~g\. 
As a result of this symmetry, the scattering state 1^'^+^) 
in Eq. © satisfies = exp{ikyd)\^^+^}. There- 

fore, we can conclude that coefficients (3j takes the form 



13 J = ^e*''"'*'. 



(10) 



Second, substituting Eq. dTU]) into Eqs. ([51 H]), we 
obtain the expression for the j-independent coefBcient /3: 



13 



1 



27r e- _ uja 



(11) 



where the self-energy S(/c) is given by 



d-l 



Here the function is defined as 



(12) 



dqx- 



1 



^0^ 



with 



Pl{ky) = {ky +TT + 



2Tr\l\ 



)mod[27r] — tt. 



(13) 



(14) 



It is pointed out that, in our system the self-energy S(fc) 
takes finite value (except for some special momentums 
which will be discussed later) and thus the renormaliza- 
tion technique is not required. That is due to the fact 
that the single-photon energy ej: has a finite upper limit 
ujc + 2^. To obtain Ea. (fT5|) we have also used the formula 



-t-oo 

E 



-i{qy-ky)dj _ 



2ti 



d-l 



-rY.^[%-Pii.ky)\^ (15) 



for Qy G [-TT, tt]. 

Furthermore, we can treat the the integral in S;(A:) 
analytically and get the result 



2dS, 



-sign(AO l^^r 2d^V^4f~T 



with Ai defined as 



Ai = cos kx + cos ky — cos[pi{ky)] 



l^il <1, 

\Ai\ > 1, 
(16) 

(17) 



Substituting the resuhs into Eqs. [51171). we 

finally obtain the analytical expressions of the state 1 0) , 
the coefficient f3j, and the scattering state 



With the analytical expression of the scattering state, 
we can calculate the on-shell element t{k' <— k) of the 
T-matrix. According to the scattering theory, t{k' <— k) 
is defined as t{k' ^ k) = The straight- 

forward calculation yields 



d-l 



t(k' ^k) = ui(fc) J2^[K- PiC'y^] ' (18) 

(=0 

where k' = (fc^, ky) and the function ui{k) is defined as 



Ul(fc) = 



\n\ 



(19) 



A — 2^ (cos kx + cos ky) ~ S(fc) 
Here the photon-atom detuning A is defined as 

A^LOc-UJa. (20) 

Due to the delta functions in Eq. ([18]), the y-component 
ky of the outgoing momentum can only take d possible 
values. That is also due to the translation symmetry 
along the y-axis of our system. 

C. Single-photon scattering probability 

Using the above results of the T-matrix element, we 
can calculate the single-photon scattering probability. To 
this end, we consider the scattering of a single-photon 
wave packet on the atoms. In the scattering process, the 
incident wave packet of the photon can be expressed as 



(21) 



Here 0*^*"' {k) is single-photon wave function in the mo- 
mentum representation and satisfies J (ifc|0*^*"-'(fc)|'^ = 1. 
We further assume tp'^^'^^k) sharply peaks at a specific 
momentum fco = (fcoa:,fcoj/)- According the scattering 
theory, when the scattering process is completed, the 
single-photon state can be expressed as 

|$(°"*)) = / dk\k)\g){~g\{k\S\<f^^-^) (22) 
dk(j)^°'''^k)\k} \g} (23) 



in the interaction picture. Here the S'-matrix satisfies 
{g\{k'\S\k) \g) ^ 6{k' -k)~2Tii5{€^-€p)t{k' ^ k). (24) 



Substituting Eqs. ([HIMl) into Eqs. ([221 [23]) , it is easy to 
find that after the scattering process, the incident wave 
packet splits into {2d — 1) different ones (see Appendix). 
Namely, the out-put wave function in Eq. (|23p is given 

by 
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(25) 



Here the ^-th wave packet (j)f^^\k) sharply peaks at a 
momentum ki — [kix,pi {k^y)] which satisfies e^^ — 
and sign(fci2;) = sign(Z) (see Appendix). 

Then the probability for the photon being scattered 
to the l-th {I ^ 0) outgoing momentum ki is P; — 

f dk\(j)[°^*\k)\'^ . As shown in the Appendix, Pi can be 
expressed as 
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4^2 \smkox sinfcia; 



(26) 



Therefore, for an incident photon with central momen- 
tum fcg, the scattering probability i?i(fco) is 



d-1 



l=~{d~l) 



(27) 



1=1 



where we have used the fact Pi — P_ 



D. Behavior of the scattering probabihty 

Now we discuss the behavior of the scattering probabil- 
ity i?i(fc) for a incident photon with central momentum 
k. According to Eqs. (gTl ^ and Eq. HI]), it is clear 
that for fixed value of fc, Ri{k) is a Lorentz function of the 
photon-atom detuning A, and takes the maximum value 
under the condition A = 2^ (cosfc^^ -I- cosfcj^) -I- Re[I](fc)]. 
In Fig. [2] we illustrate Ri{k) with different periods 
of atoms. Here we choose the incident momentum 
{kx,ky) = (7r/4,7r/2). The Lorentz-shape of the i?i(fc) 
is clearly shown. 

The single-peak behavior of i?i(fc) can be explained 
by the following simple picture. Owing to the periodical 
structure of our system in the y-direction, in the single- 
photon scattering process the incident state |fc) of the 
photon is coupled to the atomic spin-wave state 



\Sky) 



i 



(28) 



As a result of this coupling, the effective energy of state 
\Sky) is shifted from the bare value uJa to uja + Re[I](fc)]. 
Thus, Ri{k) takes the maximum value when the energy of 
the incident photon is resonance with the shifted energy 
U!a + Re[E(fc)] of the spin- wave state, or the condition 
eg = + Re[E(A;)] is satisfied. It is pointed out that, the 

energy shift Re[E(fc)] is related to the density of the 
atoms, and can be considered as the collective Lamb shift. 



Figure 2: (Color online) The single-photon scattering prob- 
ability Ri (k) as a function of photon-atom detuning A with 
d — 1 (purple solid line), 2 (green dash line) and 3 (blue 
dash-dotted line). Here we choose k — {kx,ky) — (7r/8,7r/4) 
and fl — 5, ^ = 1. The scattering probability Ri{k) has a 
Lorentz-shape profiles with Lamb shifts Re[E(fc)] which are 
related to d. 



The cooperative effect of the nuclear ensemble is reflected 
by the atomic-density-dependence of this collective Lamb 
shift and the width Im[E(fc)] of the peak of Ri{k). 

Ill Fig. [3] (a, b) we plot the scattering probability Ri{k) 
as a function of cosfc^; for fixed values of A and ky. It 
is shown that, for some certain momentum kx, we have 
Ri{k) = 0. Namely, in these cases the photon is not scat- 
tered by the atoms. This phenomenon can be explained 
as follows. Due to the periodical structure of our system, 
for the incident photon with momentum fc, only the pho- 
ton states |fc;) = \kix,pi (ky)), which satisfies e^: = e^, 
are involved in the scattering process. Eq. (fT3|) shows 
that, the coupling between the atoms and the photon in 
the state \ki) with a special I contributes the self-energy 
E;(fc). With straightforward calculation, we can re- write 
E,(fc) defined in Eq. 1^ as 



2 

27rd 



dx- 



p{x) 



'^P^'^^ + 2^ 



X + iO- 

2 pbi 
P 



rf,^, (29) 



ai 



where P means the principle integral, a; = £(o,pi(fcy)) 
and hi ~ ^{TT,p,{ky)) are the lower and upper bounds of 
the energy ej:^ with respect to the state |fc;), respec- 
tively. In Eq. we also have p{x) — 2dp{x)/dx with 
p{x) — arccos[x — -I- 2<^ cos\pi{ky)\\. Usually, the self- 
energy is convergent for a system with bounded energy 
spectrum. However, in the present problem, when the en- 
ergy of the incident photon is in the bound of the energies 
of the states |fc;) (i.e., the condition e^: = a/ or — bi 

is satisfied), the self-energy I]/(fc) diverges according to 
Eq. (HH). Furthermore, Eqs. dill US]) and Eq. ^ yield 
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Figure 3: (Color online) (a) and (b): The single-photon scattering probability Ri{k) as a function of coskx with ky = 7r/4 (a) 
and 7r/3 (b). (c) and (d): the real part (red solid line) and imaginary part (black dashed line) of self-energy E(fc) as functions 
of kx with ky — 7r/4 (c) and vr/S (d). Here, we also choose f2 = 5, ^ = 1, A = and d = 3. 



i?i(fc) = when the self-energy diverges. 

The above observation is also verified by Eq. 
which shows that Ei(fc) = oo when = 1. In addi- 
tion, we illustrate the self-energy S(fc) as a function of 
kx in Fig. [3] (c, d). It is clearly shown that S(fc) = oo 
when Ri{k) = 0. 



III. SINGLE-PHOTON SCATTERING WITH 
TWO LAYERS OF ATOMS 

In the above section we have studied the single-photon 
scattering on one layer of atoms in a 2D cavity array. We 
show that the single-photon scattering probability takes 
maximum value when the incident photon is resonant 
with the shifted atomic energy Ua + Re[E(fc)]. In this 
section, we consider the single-photon scattering in the 
2D cavity array with two layers of atoms located in the 
cavities at {xi,dj) and {x2, dj) with j = 0, ±1, ±2, ... (Fig. 
|4|). We will show that as a function of A, the scattering 
probability has two peaks rather than a single one. The 
double-peak behavior results from the photon-induced ef- 
fective coupling between atoms in different layers, and 
can be considered as an EIT-like phenomenon. 

In the presence of two layers of atoms, the atom-photon 




Figure 4: (Color online) The single-photon scattering by two 
layers of atoms in a 2D CCA. In this case the atoms are 
confined in two layers at a; = a;i and x = X2, with the same 
periodic d. 

interaction reads 

oo 

^ - E E ^^sH^.M) ' (5l + h.c), (30) 

s— 1,2 j — — oo 

with |g(e))j'''' the ground (excited) state of the j-th atom 
in the s-th layer. The scattering state 1^'^+^) can be 
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written as 

|^(+)\ 



E E /3f'|vac)|etO, (31) 

s— 1,2 J — — oo 



with \g) = Y[j \9)f'' collective ground state 



of all atoms, and |ej'*'') (s = 1,2) defined as je^*'') = 

riiT^j Is);"'' rij denotes the state in which 

only the j-th atom in the s-th layer is excited. More- 
over, the translation invariance along the y-axis leads 
to the result = (3^'^ exp[i{kxXs + kydj)]. Substitut- 
ing this result into the Lippmann-Schwinger equation, we 
find that the j-independent coefficients /S*-^'^^ have sim- 
ilar expressions with the parameter /? in Eq. (|lip . and 
can be written as 



27r 



- CJa - S(fc) 



^2 



(32) 



Here the self-energy S(fc) is now a 2 x 2 matrix 
S(fc) 

with elements (i, j — 1, 2) given by 



S21 ^22 



(33) 



1=0 



E 



exp[-z(fca; - qx){xi - Xj)] 



ffe p, (fey)] +*0^ 



(34) 

It is clear that Eq. ([5^ can be solved straightforwardly 
and we have 



27r(I]+-S- 

for s = 1,2. Here = ujj: — uja and is defined as 



, (35) 



J, = S 



17 



(3-s),(3-s) 



{3-s) . 

17,, 



,(3- 



(36) 



In Eq. S±(A;) is the eigenvalue of matrix S(fc) and 

takes the form 



1 



Ell -f E22 ± Y (Sii — S22) +4E12S21 



(37) 

With these results, we can derive the expressions of the 
scattering state and the on-shell element t{k' •<— 

k) of the T-matrix is 

t{k' ^k)= uu{k) Y.S[k'y- PI {ky)] . (38) 

Now the coefficient function uii(A:) is given by 

uu{k) = i y e-'('=^-'=-)"=17:/3(^). (39) 

s=l,2 



Similar as the above section, for an incident photon 
with central momentum fc, the scattering probability 
Ru{k) can be expressed in terms of the T- matrix ele- 
ment, and the straightforward calculation gives 



Ruik) = 2j2 



uu{k) 



4^2 \smkx sinfc;^;! ' 



(40) 



with kix determined by the equation = • 

Now we investigate the behavior of the scattering prob- 
ability -Rii(fc) with respect to the bare detuning A. With 
Eq. and Eq. ^5^, we find that -Rii(A:) takes local 



maximum values when the condition Re[Aj: — 'E±(k)] — 
is satisfied. Namely, unlike the single-peak behavior 
shown in Eq. (|lip and Fig. ^ for the single-layer case, 
in our current system Ru{k) has two peaks around the 
positions 



A± = Re[E,±] + 2^ (cos k^ + cos ky 



(41) 



Furthermore, in the two-layer case, the cooperative effect 
of the nuclear ensemble is reflected by the atomic density 
dependence of the Lamb shifts Re[I]±(fc)], the peak 
widths Im[S±(A;)] and the distance Re[E+(A:) — E_(fc)] 
between the two peaks. 

This observation is verified by our exact numerical cal- 
culation for Eq. (gni)- In Fig- ®, we illustrate Ru{k) 
for the one- and two-layer cases. The double-peak be- 
havior of i?ii(fc) at A = A± for the two-layer case is 
clearly shown in Fig. ([S^-lSl:). It is pointed out that, such 
a behavior essentially has the same physical mechanism 
with that of the atomic susceptibility in an EIT system. 
In typical EIT system, two internal states of the A-type 
atom are dressed with the control laser beam, and form 
two non-degenerate dressed states. Then the atomic sus- 
ceptibility takes local maximum value when the incident 
photon is resonance with one of these two states. In our 
problem, the incident state |fc) of the photon is coupled 
to two quasi spin- wave states 



(42) 



with respect to the excitations of the atoms in the 1st 
and 2nd layers, respectively. These two states have the 
same bare energy uja- Nevertheless, due to the atom- 
photon coupling, isj^'') and |5'^^'*) are effectively coupled 
with each other via the non-diagonal elements S12 and 
S21 of the self-energy matrix S(fc), and form two dressed 
states [Fig. These two dressed states have different 
effective energies Wa + The single-photon scattering 
probability i?ii(fc) takes local maximum values when the 
incident photon is resonant with one of these two dressed 
states. Therefore, the i?ii(fc) presents an EIT-like behav- 
ior which is illustrated by its two peaks at A = A±. 

In the end, we emphasis that the above analysis is a 
qualitative one which provides a physical picture for the 
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Figure 5: (Color online) The single-photon scattering probability -Rii(fe) for the two-layer case (solid blue line). The peak 
behavior of Rii{k) at A = A± is illustrated for X2 ^ xi — 15, {kx,ky) — (tt/S, 7r/4) (a), X2 — xi = 14, (kx,ky) = (tt/S, 7r/4) 
(b), X2 — x-i = 12, {kx,ky) — (tt/S, 7r/3) (c) and X2 ~ xi = 8, {kx,ky) — (tt/S, 7r/3) (c). Here we choose ^li — fl2 = 5, ^ — 1 
and d = 3. As a comparison, we also plot the scattering probability Ri (k) for the single-layer cases with the same parameters 
(dashed red line). 



(a) (b) 




Figure 6: (a): In the two-layer case, the two quasi spin- wave 
states 1 5^^') and IS^^'), which have the same bare energy uja, 
are coupled to the continuous spectrum of the photonic states. 
Due to this atom-photon coupling, and |S'^'^') are effec- 

tively coupled with each other by the non-diagonal terms E12 
and E21 of the self-energy matrix E(fc). (b): As a result 
of this effective coupling, IS^}^) and IS^^') form two dressed 
states with different effective energies uja + S± . 

EIT-like behavior of the single-photon scattering proba- 
bility i?ii(fc). On the other hand, due to the complicated 
expressions of the complex self-energy matrix S(/c), in 
different parameter regions, the widths and heights of 
the two peaks of Ru{k) can be quite different, as shown 
in Fig. ([S]). Moreover, in some special cases the peak- 
widths could be larger than the inter-peak distance, and 



thus the two peaks emerge into a single one Fig. ([SJl). 



IV. CONCLUSION 

In this paper, we study the single-photon scattering 
by one or two layers of identical atoms in a 2D CCA. 
We obtain the exact analytical expression for the single- 
photon scattering probability, and explore the relevant 
cooperative effects. In the case of one atomic layer, a 
Lorentz-type peak of the scattering probability implies 
the collective Lamb shift, while in the two-layer case an 
EIT-like double-peak behavior appears, and the cooper- 
ative effect of the ensemble is reflected by the atomic- 
density-dependence of the position, width and distance 
of the two peaks. Our approach fully characterize the 
quantum nature of the photon and the internal structure 
of the nucleus, and will be helpful for the study of the 
2D photonic quantum devices and X-ray quantum optics 
from the perspective of quantum scattering. 
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Appendix: single-photon outgoing wave packets and 
probability Pi 

In this appendix derive the outgoing wave function of 
a single-photon scattered by the atoms, and prove Eq. 
(HSl) for the probability Pi. Substituting Eqs. (HH ^ 
into Eqs. it is easy to find that after the scatter- 

ing process, the incident wave packet splits into {2d — 1) 
different ones. Namely, the out-going wave function in 
Eq. (HSl) is given by 



Here the function (j)[°"*\k) is defined as 

~ 2^|sinfc,| ' 



(43) 



for / = and 



^(out)^^^^ u[fi{m^"^)[fi{k)] 



2C sin/ia;(fc) 
with the function /;(fc) = fix{k), fiy{k) 



^sign(A;^),sign(/) 



(44) 



(45) 



defined as 



//. (fc) 



= arccos -j cos k^ + cos ky — cos 



//.(fc)]}, 

(46) 

fiv(f) = (fc, + ^+^)mod[27r]-7r, (47) 



and the Kronecker symbol 5i^j satisfies 5i^j = 1 for i =^ 
j and 5i,j = Q ioi i ^ j. In Eqs. (gi]) and (gS]), we 



have u{k) = ui{k) for the case with one layer of atoms 
and u{k) = uii(k) for the two-layer case. Then it is 
easy to see that, when the incident wave packet 0'-™^(A:) 
sharply peaks at a specific momentum /cq = {kox,koy), 



the out-put wave packet 



Aout) 



{k) sharply peaks at the 



momentum fc/ = [kix , pi {k^y)] which satisfies e^^ — ej:^ 
and sign (kix) = sign (l). 

Now we calculate the probability Pi — f \4>\°^^\k)\'^ dk 
for the cases with 1^0. Apparently, we have 



P 



u[/7 (fc) [/7 




2^ 


sin fix (jtj 





^sign(fc^ ),sign(/) 



dk (48) 



We define gix — fix{k) and giy — fiy{k). Then using the 
fact that (/)*^"') (fc) sharply peaks at a specific momentum 
ko and the relations 



dkx 



sin- gix 



sin kx 



dgix , dky — dkiy , 



(49) 



we have 



P 



u[fi{kM(^^)[fi{k)] 



2i 



sin fix{k) 





u{ko) 


2 


4e\s 


m kox 
u{ko) 


3in kix 1 

2 



4^2 Isinfcoa; sin kix \ ' 
Then we have proved Eq. (P^ . 



dk 

2 

(f'^'"H9ix,9iy) dgixdgiy 
(50) 
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